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ABSTRACT 

We consider how upcoming photometric large scale structure surveys can be opti- 
mized to measure the properties of dark energy and possible cosmic scale modifications 
to General Relativity in light of realistic astrophysical and instrumental systematic 
uncertainities. In particular we include flexible descriptions of intrinsic alignments, 
galaxy bias and photometric rcdshift uncertainties in a Fisher Matrix analysis of 
shear, position and position-shear correlations, including complementary cosmolog- 
ical constraints from the CMB. We study the impact of survey tradeoffs in depth 
versus breadth, and redshift quality. We parameterise the results in terms of the Dark 
Energy Task Force figure of merit, and deviations from General Relativity through an 
analagous Modified Gravity figure of merit. We find that intrinsic alignments weaken 
the dependence of figure of merit on area and that, for a fixed observing time, halving 
the area of a Stage IV reduces the figure of merit by 20% when IAs are not included 
and by only 10% when IAs arc included. While reducing photometric redshift scatter 
improves constraining power, the dependence is shallow. The variation in constrain- 
ing power is stronger once IAs are included and is slightly more pronounced for MG 
constraints than for DE. The inclusion of intrinsic alignments and galaxy position 
information reduces the required prior on photometric redshift accuracy by an order 
of magnitude for both the fiducial Stage III and IV surveys, equivalent to a factor 
of 100 reduction in the number of spectroscopic galaxies required to calibrate the 
photometric sample. 

Key words: cosmology: gravitational lensing: weak - dark energy - equation of state 
- cosmological parameters - large-scale structure of Universe 
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1 INTRODUCTION 

The dawn of "precision cosmology" was heralded by the re- 
sults of surveys which, for the first time, produced data of 
sufficient quantity and quality that our cosmological probes 
could begin to accurately measure some of the fundamen- 
tal properties of the Universe. The outcome has been the 
ACDM concordance cosmology, a description of the Universe 
compatible with the joint constraints from many sources of 
cosmological information. The next decades will see a step- 
change in our ability to measure cosmological parameters as 
new surveys produce orders of magnitude more data than 
has previously been available. This will allow us to test the 
standard cosmological model as never before. 

The standard ACDM model describes a Universe made 
up of ~75% Dark Energy, a smooth, negative pressure fluid, 
~20% Dark Matter, collisionless massive particles which in- 
teract solely via gravity, and just ~5% bary ons making up 
the potentially visible mass of the Universe (|Dunklev et alj 



l2009t ). The standard model also assumes gravity to be de- 
scribed by Einstein's General Relativity (GR). 

A Universe governed by GR and populated by standard 
gravitating matter cannot explain the observed acceleration. 
The dark energy fluid was proposed to solve this paradox, 
reviving the idea of a cosmological constant which Einstein 
had originally included in his GR field equations and sub- 
sequently discarded. Much effort has been devoted to char- 
acterising the nature of dark energy- particularly attempts 
to discriminate between a pure cosmological constant and 
a dynamic scalar potential with time-varying equation of 
state w{a) (|Albrecht et al.|[2006t ). In all its forms, dark en- 
ergy poses problems of fine-tuning for which we have, as yet, 
no physical motivation. 

Rather than invoke dark energy, an alternative expla- 
nation for cosmic acceleration has been proposed - that our 
standard theory of gravity is incomplete and that a cor- 
rect theory of gravity would explain cosmic acceleration 
at late times and large scales in a universe populated by 
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matter with positive pressure. There are a large number 
of theoretically motiv ated modified theories of gravity, see 
iJain fc Khour"vl d20ich for a review. Testing the theory of 
gravity on cosmic scales will be one of the most interesting 
opportunities afforded by upcoming survey data. 

Weak gravitational lensing (WL) is a particularly use- 
ful probe of gravity because it is sensitive to <f> + ip, the sum 
of the metric potentials. Through the modified growth of 
structure and geometry via the lensing integral cosmic shear 
can constrain both the ratio of metric potentials and mod- 
ifications to the Poisson equation. In contrast probes such 
as galaxy redshift surveys and galaxy peculiar velocities de- 
pend only on the Newtonian potential, V- WL constraints 
also have the benefit that they probe the dark matter distri- 
bution directly, avoiding the impact of galaxy biasing. Joint 
constraints from combining multiple probes can help break 
degeneracies between free parameters, producing the tight- 
est constraints. 

Cosmic shear is the name given to weak gravitational 
lensing in random pat ches of the s ky and was first de- 
tected a decade ago ([K aiser ct all l200d: IWittman et all 



l200d ; Ivan Waerbeke et all l200fj| ; iBacon et alj 120001 ). The 
latest constraints on cosmology come from the Hubble 
Spac e Telescope COSM OS survey ijSchrabback et all 
120091 ; iMassev et all 120071 ), and the Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS) |Fu et alj|2008h. In ad - 
dition the 100 square degree survey ( Benjamin et al.ll2007l) 
comb in es data from several s malle r s urveys dHoekstra et al.l 
20021; iHetterscheidt etafl 120071 ; iLe Fevre et all 12004 



Hoekstra et al.ll2002T ). 



Cosmic shear has been identified as the method with 
the most potential to uncover the natur e of dark energy 
|Albrecht et all 120061 ; iPeacock etahl 120061 ) and therefore a 
number of surveys are planned with a major cosmic shear 
component. Upcoming "Stage III" projects include Kilo- 
degree Survey (KIDS) on the Very Large Telescope (VLT) 
Survey Telescope (VST), the Panoramic Survey Telescope 
and Rapid Response System (Pan-STARRS) project, the 
Subaru Meaurement of Images and Redshifts (HSC) sur- 
vey using HyperSUPRIMECam, the Dark Energy Survey 
(DES) on the Blanco Telescope. More ambitious "Stage IV" 
imaging projects are the Large Synoptic Survey Telescope 
(LSST) ground-based project and in space the proposed Eu- 
ropean Space Agency mission Euclid and the NASA pro- 
posed Wide-Field Infrared Survey Telescope (WFIRST). 

Greater accuracy demands better treatment of system- 
atic effects. One of the main systematics in cosmic shear 
studies is the contamination by the intrinsic alignment (IA) 
of galaxy shapes. This systematic manifests in two main 
ways, in the first physically close galaxies have preferen- 
tially aligned ellipticities, this is known as the Intrinsic- 
Intrinsic (II) correlation. In the second galaxies close on 
the sky by separated radially are anti-correlated as the 
foreground gravitational potential shapes the closer galaxy 
while also gravitationally lensing the background galaxy, 
this is called the Gravitational-Intrinsic (GI) correlation. 
The former adds to the cosmic shear signal we wish to 
mea sure, the latter subtracts from it. In an ear lier pa- 
per (|l. Laszlo. R. Bean. D. Kirk and S. Bridlell201ll ). herein 
"LBKB", we established how, for a fixed survey specifica- 
tion, the inclusion of a realistic IA model and the inher- 
ent uncertainties in the model significantly degrade the con- 



straining power of a cosmic shear survey for dark energy 
and modified gravity c osmologica l mod e ls. This work deep- 
ened t he conclusions of lBernsteinl (|2009l ) ; Ijoachimi fc Bridlel 
|2009l ) and expanded them into the MG domain. 

What is encouraging, however, is that the use of 
position-position correlations, nn, and particularly position- 
shear, ne, cross-correlations can go some way towards mit- 
igating the impact of IAs. This data is already collected 
by a standard WL survey. In the presence of IAs, con- 
straints from ee + ne + nn are roughly twice as strong as 
those from ee alone. Until now, the discussion of optimal 
survey parameters has generally assumed that ee correla- 
tions alone will be included, without IAs. In this paper we 
investigate the impact of survey strategy and design on dif- 
ferent combinations of probes, with and without IAs. This 
is a continuation of the cosmi c shear survey optimisatio n 
work begun by , amo n g others, Amara fc Refregierl d2006h : 
Bridle fc I Kinel (120071) ; J joachimi fc Bridlel d200g); iMa etall 



(2ooe 



2006); Huterer et all 



In this paper we expand on the work on LBKB to fully 
understand how the observing strategy of a given survey is 
central to the type of data and the quality of results that sur- 
vey will produce. Different survey geometries with the same 
instrument will produce different statistical errors due to a 
different balance between survey characteristics, the most 
important for cosmic shear being survey area galaxy num- 
ber density on the sky and median redshift of the galaxy 
distribution. Other properties such as required redshift ac- 
curacy are important guidelines for instrument designers 
and those prepa ring analysis pipelines and follow-up or cali- 
bration studies |Amara fc Ref rcgicr 2006]). While more area 
and increased number density will always be desirable, it 
is important to keep firm goals in mind when producing 
desiderata for future surveys which will always be limited by 
technology and finite observing time. In particular learning 
about dark energy and modified gravity may benefit from 
different survey strategies and call for the prioritisation of 
different properties. 

The paper is organised as follows: in section [2l we sum- 
marise our standard cosmology, fiducial surveys and models 
for deviations from GR, Intrinsic Alignments and bias pa- 
rameterisation, as well as our figures of merit for DE and 
MG. In section Owe investigate the sensitivity of dark en- 
ergy and modified gravity constraints to changes in survey 
specification. In section [37T1 we present the effect of varying 
survey area and the impact of finite survey time on overall 
strategy and the ability to constrain MG. The importance 
of photometric redshift accuracy is addressed in sections [32] 
and 13.31 Conclusions are made in section 2] 



2 COSMOLOGICAL SET-UP 

The paper deals with the impact of various aspects of sur- 
vey strategy and redshift quality on the power of cosmic 
shear and galaxy position information to constrain dark en- 
ergy and deviations from General Relativity. In 12. f I we sum- 
marise our cosmic shear formalism, I2.2l extends it to include 
intrinsic alignments and 12.31 adds galaxy position auto- and 
cross-correlations and introduces a coherent biasing formal- 
ism for IAs and galaxy bias. 1 2. 4 1 re views the MG parameter- 
isation we use and how it enters our angular power spectrum 
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integrals. In I2.5l we summarise some basic cosmology we use 
throughout the paper and 12.61 describes the fiducial survey 
parameters we use in the following sections. 



The total observed lensing signal is then the sum of the 
cosmic shear and the IA terms 



(4) 



2.1 Cosmic Shear 

Cosmic shear is the shape-distortion induced in the image 
of a distant galaxy due to the bending of its light by gravity 
as it passes massive structure in the Universe. If we assume 
General Relativity holds then we can define the cosmic shear 
angular power spectrum under the Limber approximation as 

§w i (x)w r i(x)i , «(M) 



cg°(0 



V 



(1) 



where Ps$(k,z) is the three dimensional matter power spec- 
trum, x is the comoving distance in units of fe _1 Mpc and 
W(x) i s the lensing efficiency function 



Wi(x) = ^-p{z)a\z) X 



2.2 Intrinsic Alignments 



(2) 



The intrinsic alignment (IA) of galaxy ellipticities is a prime 
contaminant to the measured cosmic shear signal. A naive 
approach to cosmic shear assumes that galaxy's intrinsic el- 
lipticities are randomly distributed on the sky so, when we 
average over observed ellipticity in a small patch, intrinsic 
ellipticity cancels and we are left with the induced shear. 

Unfortunately, this assumption is invalid because 
galaxy ellipticities are aligned due to two effects arising from 
the same physical intrinsic alignment origin. Physically close 
galaxies tend to align with the local gravitational tidal field 
and so are positively correlated with each other, this is the 
Intrinsic-Intrinsic (II) alignment. Background galaxies can 
have their light lensed by foreground gravitational fields 
which align the intrinsic ellipticity of foreground galaxies. 
This induces an anti-correlation and is the gravitational- 
intrinsic (GI) alignment. 

We follow the procedure of LBKB and implement IAs 
using an updated ver sion of the corrected LA model of 
iHirata fc Seliakl (|2004h . incorporating their correction of er- 
ratum 2010 and assuming that all I A physics occurs at a 
high-z epoch of galaxy formation hence the IA signal de- 
pends only on the linear matter power spectrum rather 
than the subsequent nonlinear evolution. The rest of the 
LA model enters as a factor of 



-Cipm{z = 0) 



(3) 



where p m (z = 0) is the matter density today and C\ = 
5 x lO~ 14 (/i 2 M Mpc" 3 ) _1 i s the amplitude of th e IA term, 
normalised to redshift zero (|Bridle fc Kinell2007D . This fac- 
tor appears in the IA projected angular power spectra, 
linearly in the Gravitational-Intrinsic (GI) correlation and 
squared for the Intrinsic-Intrinsic (II). The window function 
associated with IAs is the galaxy redshift distribution, n(z), 
see table [1] for the full angular power spectra equations. 

In section l2"31 we expand our notation to allow nuisance 
parameters to parameterise our knowledge of IAs, galaxy 
bias and their cross-correlations. In this more general case 
bi given in eqn.[3]becomes the fiducial value of a variable IA 
term which appears quadratically in the II power spectrum 
and linearly in the GI and gl power spectra. 



2.3 Galaxy Position Data 

A cosmic shear survey contains galaxy position information 
(angular position on the sky and redshift) as well as mea- 
surements of galaxy shear. [joachimi fe B ridle (200jJ) provide 
a formalism for including this additional information in cos- 
mological parameter constraints and show how this extra 
information can serve to partially mitigate the impact o f 
IAs. This approa ch follows from the work of IZhanel (|2008l ); 
iHu fc Jain! (|2004 ) 

The extra observables we use are galaxy position- 
position power spectra and the position-shear cross-spectra, 
defined analogously to the shear-shear power spectrum: 



QJM _ Qgg 



C< 



9 ft 



(5) 
(6) 



We now assign nuisance parameters to each component C; 
in a self-consistent way as explained in LBKB. The four 
bias function are b g , bi, r g and rj. b g and bi model bias 
from galaxy position and IAs amplitudes, and their cross- 
correlations are modelled by r a and 77. This u nified ap- 
proach was first introduced by iBernstein] |2009l ). Table [T] 
summarises the full power spectra equations, consistent with 
the definitions in LBKB eqns. 37, 38, 43 and 44. 

We let each bias parameter vary in amplitude and as a 
function of scale and redshift using a Nk x N z grid of free pa- 
rameters interpolated over fc-,z-space, i.e. each nuisance fac- 
tor, bx = A x Qx(k, z) is the product of a variable constant 
amplitude parameter, Ax, and a variable grid, Qx(k,z) in 
k and z. Throughout this paper the grid size is set to the 
fiducial value of Nk = N z — 5, which means marginalisation 
over 104 nuisance parameters when the full ee+ne+nn probe 
combination is considered in the presence of IAs. The nui- 
sance parameters multiply linearly into the angular power 
spectra integrands and each power spectrum depends on a 
subset of the bias parameters as follows: 



(7) 



As mentioned in the previous section, the fiducial value of 
bi is given by eqn.(3] b g has fiducial value 1 while r g and n 
vary around the fiducial value 0.9 to avoid more than perfect 
cross-correlation. 

We have ignore d the effect of lensing mag nification and 
followed LBKB and I Joachimi fc Bridle! (|2009h in applying a 
cut on multipole I to any redshift bin combination ij which 
includes galaxy position information (i.e. nn or ne) accord- 
ing to i max (i) = This aims to account 

for uncertainties in the galaxy bias model at small scales. 
For ne bin pairs, the cut is made on the galaxy, i.e. n, bin. 
For position-position, nn, pairs t here is a choice of bin on 
which to apply the cut. We follow [joachimi fc Bridle] (|2009T 1 
in making the optimistic choice and cutting on the higher 
redshift bin. Shear-shear, ee, pairs do not depend on galaxy 
bias and are therefore used up to the full range in I. 
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Table 1. Summary of the projected angular power spectra considered in this work. 
Correlation 2D PS 



Shear 

Intrinsic-shear 
Intrinsic 

Galaxy clustering 
Clustering-shear 

Clustering-intrinsic 
Galaxy ellipticity (observable) 
Galaxy number density (observable) 
Number density-ellipticity (observable) 



cg G (Q 
eg 1 ® 



4(0 



= /^.(x)^(x)[QW^] Pss(k,z) 



^ MMW (M) 

/ #m(x)%(x)^(fe, z)Q2(z f )R2(z f )P ss (k, z) 
C° J 9 (l) = f$n i (x)n 3 ( X )b 2 g (k,z)P M (k,z) 

Cf(l) = f ^n l ( X )W j ( X )Q(z) 1 -±^b g (k,z)r a (k,z)Pss(k,z) 



/ ^n r { X )n 3 ( X )Q(z f )R(z f )b g (k, z)bj(k, z)r g (k, z)rj(k, ^ ^Wg g 



CP. 1 



C: 



go 



2.4 Modified Gravity 

There are a large number of modifications or extensions 
of Einstein's general theory of relativity on cosmic scales 
which come under the general heading of Modified Grav- 
ity theories. These can be motivated by the presence of 
extra dimensions, as in DGP, or extra degrees of freedom 
comp ared to the GR a ction equation, as i n f(R), TeVeS 
etc. (iDvali et all l200d ; I Carroll et all l200rj ; ISkordij 120091 ; 
J am fc Khourvll2010l ). As in LBKB, we do not assume a 
particular modified theory of gravity but rather concentrate 
on "trigger parameters" whose deviation from their GR val- 
ues would indicate the presence of some physics beyond that 
in the standard GR picture. 

In the conformal newtonian gauge the metric for a flat 
FRW spacetime is written 



Fiducial Survey Paramters 



ds 



-a(-r) 2 [1 + 2tp(x, t)] dr 2 + a(r) 2 [1 - 2<f>(x, i)] dx 2 



(8) 

where if) and <j> are the scalar potentials which describe per- 
turbations to the time- and space-parts of the metric respec- 
tively. We parameterise deviations from GR through two 
parameters, Q and R. One alters the way the Newtonian 
potential responds to mass via the Poisson equation, 

k 2 ip(x,t) = -4-KGQpa 2 8 (9) 

and the other modifies the ratio of the metric potentials 

ip(x,t) = Rtf>(x,t) (10) 

Q and R are assumed to be scale- independent and vary with 
redshift as 



Q = (Qo - l)a s 
R=(Ro- l)a s 



(11) 
(12) 



where a is the scale- factor and Qo,Ro are the free parame- 
ters we vary for MG, and s — 3. We are interested in MG 
theories which explain the acceleratin expansion of the Uni- 
verse observed at late times, s = 3 allows any modification 
to "turn-on" and late times and avoids violating early Uni- 
verse constraints from the Cosmic Microwave Background 
(CMB) and Big Bang Nucleosynthesis (BBN). 

The modified gravity parameters enter the projected 
angular power spectra in different combinations, as shown 
in Table [T] due to their different dependencies on the mater 
density field. As well as modifying the angular power spec- 
tra, Q and R enter the angular power spectrum integrals 
via their response to the metric potentials through the 



Parameter 


Stage III 


Stage IV 


Area 


5,000deg 2 


20,000deg 2 


n g 


10 


35 




0.23 


0.35 


N z 


5 


10 


S z 


0.07 


0.05 


feat 








A z 


1 


1 


Zf) 


0.8/1.412 


0.9/1.412 


a 


2 


2 


P 


1.5 


1.5 



Table 2 . Fiducial survey parameter values for the DES-like (stage 
III) and Euclid-like (stage IV) surveys. 



linear growth function. The potential for weak lensing to 
co nstrain deviations from GR h as previously bee n note d 
in iBean fc Tangmatithaml <|201Ch : iLaszlo fc Beanl (|2008h : 
iBevnon et al.l (|2009l ). 

The full projected angular power spectra, including IAs, 
MG and bias parameters, is summarised in Table [1] Note 
also that deviations from GR will also manifest in the growth 
function which produces the matter power spectra in Ta- 
ble [T] We include modifications to growth via a ratio of 
MG/GR power spectra calculated using the modified ver- 
sion of CAMB used for LBKB. For ease, LBKB also provide 
a fitting function to compute this ratio over a range of Qo, 
Ro and s values. 



2.5 General Cosmology 

Throughout this paper we assume a flat ACDM cosmological 
model with fiducial parameter values equal to the WMAP7 
best-fit values: Q m = 0.262, tt b = 0.0445, w = -1, w a = 0, 
a 8 = 0.802, h = 0.714, n 3 = 0.969, Q„ = 0. Where fl m ,n b 
and Q v are the dimensionless matter, baryon and neutrino 
densities respectively, wp and w a are the dark energy equa- 
tion of state parameters l|Albrecht et al.ir2006h .o« is the nor- 
malisation of the linear matter power spectrum, h is the di- 
mensionless Hubble parameter today and n s is the power law 
of the primordial power spectrum. The lin ear matter power 
spectr um is given by the fitting formul a of Eisenstein fc Hul 
(|1998T I. with nonlinear corrections from lSmith et al I (|2003l ). 
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When we treat MG we assume that the background 
expansion of the Universe is consistent with our fiducial dark 
energy model and deviations due to non-GR physics enter 
through the growth of structure from initial perturbations. 
Unless otherwise specified all plots include CMB priors from 
a Planck-type survey as described in LBKB. We divide our 
survey galaxy redshift distributions into redshift slices with 
equal galaxy number density to allow redshift tomography 
where auto- and cross-correlations of these redshift bins are 
considered, allowing us to measure redshift evolution of the 
observables. 

We calculate constraints on cosmological parameters us- 
ing the Fisher Matrix formalism 



JV rf NP 



d%Jf) 



Cov m * (0 



dp v 



(13) 



where m, n label tomographic redshift bins, 2> m (l) is the 
data vector, in our case some combination of the angular 
power spectra C", C™ E and C ; " n . Cov mn (0 is the covari- 
ance matrix, defined as in ljoachimi fc Bridle! (|2009T ) . p M are 
the parameters varied and Nd is the number of indepen- 
dent combinations of tomographic bins. Unless otherwise 
stated the Fisher matrix varies the cosmological parameters: 
p = {£lm,wo,w a ,h,a 8 ,Qb,n a }. 

Using the Fisher matrix formalism the lower limit on 
the minimum variance bound on the error of a parameter p M 
marginalised over all other parameters of interest, is given 
by c(Pii) = v {F^ 1 )^^. We quote results in terms of the 
DETF FoM for Dark Energy 

1 



FoMde 



\/det [(F G r) W0 ,u 



(14) 



Here (Fq R ) W0:Wii is the 2x2 submatrix of the inverted Fisher 
matrix for cosmological and nuisance parameters excluding 
the MG parameters are excluded (these are fixed at their 
GR values). The FoM is proportional to the inverse of the 
area of the constraint contour in wq — w a space. 
By analogy we define a MG FoM as 



FoMmg 



det [(F 1 )q ,q (i+r /2) 



(15) 



where we have inverted the full fisher matrix including the 
modified gravity parameters and dark energy parameters. 

2.6 General Survey Parameters 

LBKB investigated the degeneracy between IAs and MG pa- 
rameters for a fixed survey specification. It was found that 
the inclusion of a realistic IA model reduced the FoM (dark 
energy or MG) for a typical stage-IV cosmic shear survey 
by ~70%. The effect could be substantially mitigated by 
the inclusion of galaxy position data and galaxy-shear cross- 
correlations. The constraining power of ee + ne + nn, includ- 
ing IAs, is roughly double that of ee, including IAs, alone. 
We showed that our results were robust to different num- 
bers of free nuisance parameters, accounting for flexibility 
in the galaxy clustering and IA bias models. The results 
were broadly similar for attempts to constrain dark energy 
equation of state parameters and deviations from GR. 

In this paper we vary certain survey parameters around 



a fiducial cosmic shear survey set-up, corresponding to Stage 
IV project. This survey has an area = 20000deg 2 , number 
density of galaxies projected onto the sky, n g — 35 arcmin~ 2 , 
and galaxy intrinsic ellipticity dispersion, a 1 — 0.35, a Gaus- 
sian photometric redshift scatter of width S z — 0.05 which 
is related to the (redshift dependent) rms photometric dis- 
persion via a z (z) = (l + z)8 z , with a fraction of catastrophic 
outliers, / cat = 0. We analyse the results using cosmic shear 
tomography with 10 tomographic redshift bins of equal num- 
ber density. 

The redshift distribution of galaxe s, n(z) is assumed to 
be given by a Smail-type distribution (|Smail et al.lll994h 



n(z) oc z a exp 



2() 



(16) 



with a = 2, P = 1.5 and zo = Zm/1-412, where the median 
redshift, z m = 0.9. 

For comparison purposes, some results are also pre- 
sented for the stage-Ill Dark Energy Survey (DES), due to 
see first light in early 2012. Table [2] summarises the specifi- 
cations for both surveys. 



3 SURVEY SPECIFICATIONS AND 
FORECAST CONSTRAINTS 

In this section we analyse the sensitivity of the prospective 
dark energy and modified gravity figures of merit on vari- 
ations in survey specification, principally focusing on their 
relationship to the survey area and depth, and the precision 
of the photometric redshift measurements when uncertain- 
ties in how intrinsic alignments are included in the modeling. 

Throughout, we quantify constraints using the dark en- 
ergy (DE) and modified gravity (MG) figures of merit in 
(|14[) and (|15l) normalised relative to a 'baseline' figure of 
merit for the shear-only auto-correlations l GG\ with IAs 
excluded, using the fiducial survey described in Table [2] 



3.1 Survey area 

Figs. [T] and [2] present the key results showing how the con- 
straining power of a survey, for dark energy equation of state 
and modified gravity parameters, varies with the survey 
area. We consider two alternative approaches, firstly a sim- 
ple option where the increased area derives from an increase 
in survey time, and secondly a more realistic option when 
survey time is fixed, and increased survey area is achieved 
by a reduction in survey depth/ limiting flux. 

Let us first consider the case where all survey parame- 
ters are kept fixed except for the survey area. We focus on 
the quantitative implications for the Stage IV survey, the 
results for a Stage III survey, presented in figure [2] are qual- 
itatively similar. 

When IAs are excluded from the analysis the improve- 
ment in FoM with increasing area can be roughly described 
by a power law dependence FoMocArea 1 with x = 0.4 and 
0.7 for CMB+shear correlations alone for the DE and MG 
FoMs. The MG FoM is slightly more sensitive to changes in 
area than t he DE FoM. The DE inde x is lower than that 
reported in lAmara fc Refregierl (|2006l ) only because of the 
inclusion of, survey area independent, CMB priors. 
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No lAs 



No lAs, Constant Survey Time 



1.5 



1 1 



C 0.5 





Area (deg' 



0.5 1 1.5 2 

Area (deg 2 ) x i 4 



With lAs 



1.5 



0) 

C 0.! 



DE ££ 

— DE £E+nn 

— DE ££+ne+nn 
- — DE nn 


MG ££ 

— MG £E+nn 

— MGf£+nE+nn 
-■-■MG nn 








0.5 
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Figure 1. Figures showing the variation with survey area of a relative Stage IV survey FoM for dark energy (DE) parameters wg,w a 
[blue lines], and modified gravity (MG) parameters QrjjQo(l + Ro)/2 [red lines], as described in (1 1 4t and Q15t . relative to a 'baseline' 
FoM. The baseline includes shear auto-correlations 'GG' alone, excluding IAs, over a 20,000 square degree survey. All FoMs contain 
priors from a Planck-like CMB survey. The figures show relative FoMs for an optimistic scenario, in which uncertainities in the IA 
model have been excluded, [top panels] and a conservative scenario, in which uncertainties in the IA model are marginalized over using 
a Nfc = N z = 5 gridded bias model [lower panels]. The baseline model (relative FoM=l) is shown as a black cross in the top panels. We 
consider the impact on the relative FoM of increasing survey area, by increasing survey time for fixed limiting magnitude [left panels] and 
fixing the survey time to trade-off survey area and depth [right panels] . Four data combinations arc considered in each panel, shear-shear 
correlations 'ee' [dotted lines], galaxy-galaxy positions 'nn' alone [dot-dash], 'ee+nn' [dashed] and when shear-position cross-correlations 
'ne' are included [full lines]. For the galaxy position correlations a TVj, = N z = 5 gridded bias model is used throughout. 



As discussed in LBKB, including IAs in the analysis de- 
teriorates both the DE and MG FoMs, with the MG FoM 
more negatively impacted. For the fiducial Stage IV survey 
the FoMs are reduced by 60 and 65 % relative to when IAs 
are excluded, consistent with LBKB. When IAs and their as- 
sociated uncertainties are marginalized over, the power law 
scaling is weaker, with x — 0.35 for DE, so that the relative 
degradation of the FoM becomes more pronounced as survey 
area increases. This reflects that the IAs are astrophysical 
systematics that are not significantly removed by increasing 
the survey area. 

Because of the uncertainties in the galaxy bias, through 
the gridded bias and cross-correlation coefficients, b g and 
r g , the cosmological information from the galaxy position 
data is massively suppressed. For the fiducial full survey area 
FoMs for CMB+nn are 22 and 1% of the CMB+pure lensing 
signal, for the scenarios with and without IAs respectively. 

If one adds in galaxy position autocorrelations to shear 
autocorrelations, 'nn + ee', we find a T2 and 36 % improve- 
ment in the DE and with and without IAs respectively, and 



19 and 40% for MG FoMs. While the cross-correlations, 
'nn+ne+ee', give only a small improvement when IAs are ex- 
cluded, they provide an 34% improvement on 'rm + ee' when 
uncertainties in the IAs are included. This results from the 
ne correlations breaking the degeneracies between the cos- 
mological and IA amplitude nuisance parameters. 

In reality most surveys have a fixed observing time. This 
inevitably leads to a trade off between total survey area and 
survey depth, which impacts both number density of lensed 
galaxies, n g , and median redshift achieved, z m . To imple- 
ment the fixed survey time option, we use a fitting formula 
for the number of galaxies per unit redshift per square ar- 
cminutc 



3z 

n(z) — E x —3 exp 
Zz n 



3/2' 



(17) 



(e.g. iBaugh fc Efstathioul l|l993M and we interpolate the 
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Figure 2. As in figure[T]but for a Stagc-III large scale structure survey with fiducial survey area of 5000 square degrees. 



numbers from Table 1 of iBlake fc Bridle! (|2005l ) using 



2 = 0.055(r n 
35400 



E 



00-' 



Tlim 



24) + 0.39 

19 



(18) 



where ru m is the limiting apparent magnitude in the SDSS r 
filter and Eo is an overall surface density in units of deg -2 . 
We use to scale the number density of galaxies with 
area relative to the fiducial values. We estimate the change 
in the limiting magnitude, A mag , with area in the usual way 



1.25logia 



A 



(19) 



As one might expect, the benefits of going to larger 
survey area are now less strong because of the lost depth, 
with the shear-only power law dependence weakening to x = 
0.3 for DE and 0.35 for MG. The baseline result for cosmic 
shear analyses is that it is still better to go to larger survey 
area d espite the lost depth. Figure 4 of|A mara fc Refregierl 
suggests that, in the absence of a CMB prior and 
IA modeling uncertainties, the FoM is roughly proportional 
to the square root of the area for fixed observing time. We 
find that with the inclusion of CMB priors, we obtain a 
similar result although the dependence on area is slightly 
weakened as in the constant depth case. When the galaxy 
position information is added the same story holds, but the 
dependence on area is slightly stronger. 

When intrinsic alignments are added the improvement 



with increased area is significantly reduced, particularly in 
the case of ee-only which reduces to a power law of x ~0.15 
for DE. A survey with an area of 5, OOOdeg 2 gathers ~ 80% 
of the DE information of a ~ 20, OOOdeg 2 survey. This is be- 
cause intrinsic alignments dominate over cosmic shear at 
very low redshift, and by marginalising out the intrinsic 
alignment nuisance parameters we effectively remove the in- 
formation at low redshift. Therefore we expect any cosmo- 
logical figure of merit to eventually reduce in the (unphysi- 
cal) limit of very large area and small depth. 

The galaxy clustering alone constraints are insensitive 
to the intrinsic alignments and therefore increase as usual. 
When using the full ee + ne + nn dataset there is more ben- 
efit to increased area even with the tradeoff of shallower 
depth, although in this case a 5000 square degree survey 
still gathers ~ 70% of the information in a ~ 20, OOOdeg 2 
survey. Note that the exact details of these conclusions will 
depend on the fiducial intrinsic alignment model and the 
number of intrinsic alignment and bias nuisance parameters 
marginalised over. Without galaxy position-shear cross cor- 
relations the constraints bear out the same trend in relative 
FoM with area, just with an overall lower magnitude. 

We repeated the exercise using the default survey pa- 
rameters for a stage III survey in Fig. [5] We vary the sur- 
vey area for fixed survey time, and cover the same range 
as in Figure [1] The survey depth for each survey time is 
correspondingly shallower, so that we have 10 galaxies per 
square arcminute for the 5000 square degree survey area. 
The left-hand plots in Fig. [2] tell a similar story to those in 
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Fig. \T\ When a constant survey time is considered there is a 
more pronounced tailing off of improvement with area above 
~ 5 — 10, 000 deg 2 than was the case for a stage IV survey. 
As the area is increased above 5000deg 2 the number den- 
sity of galaxies falls from its fiducial value of 10 arcmin~ 2 . 
For larger areas the number density becomes very low and 
shot noise rapidly increases at small scales. The loss of in- 
formation from non-linear scales has a strong impact on the 
constraining power of the survey, limiting the improvement 
with increasing area. The fact that nn alone does not dis- 
play this behaviour, displaying the same behaviour for fixed 
depth or survey time, supports our analysis- the galaxy po- 
sition information is already subject to a stringent cut on I, 
removing the nonlinear information and protecting nn from 
the impact of shot noise. 

3.2 Photometric Redshifts 

Accurate redshift information is essential for weak lensing 
tomography. As survey size increases the number of red- 
shifts becomes of order hundreds of millions to billions. The 
expense of acquiring spectroscopic redshifts for this number 
of galaxies is prohibitive. As a result the next generations of 
WL surveys will rely on photometric redshifts. 

In this paper we f o llow t he photometric redshift model 
ofH mara fc Refregieil (|200 6). including the statistical dis- 
persion of measured photometric redshifts and catastrophic 
failures in redshift measurement. The statistical dispersion 
is modelled assuming that measured photometric redshifts, 
z p hot are described by a Gaussian probability distribtion 
centred on the true redshift, Zt- The width of the Gaussian 
is the photometric redshift error which we assume scales as 
(1 + z) and is controlled by a paramter 5 Z , 

Pz h Az)= 1 (20) 

In addition to this uncertainty in the photometric PDF 
galaxy redshifts can be entirely misidentified if, for exam- 
ple, the wavelength range examined is insufficient to identify 
important features in the spectrum. A misidentified redshift 
can be said to be assigned a redshift offset by an amount 
A z from the true redshift. Affected galaxies are known as 
catastrophic outliers. 

We construct a PDF to describe these catastrophic out- 
liers, P Zcat , and combine it with the scatter PDF, writing the 
full probability distribution 

P{z phot \zt) = (1 - fcat)P* Fho M) + (fcat)P* eat (zt) (21) 

where f ca t is the fraction of our measured galaxies which 
suffer catastrophic redshift estimation failures. 

The distribution of P Zcat is bimodal, reflecting the fact 
that galaxies can either be misidentified as being higher or 
lower than their true redshifts, 

< P z ± t >= zf at = zt ± A z (22) 

where P Zaat = P - + P + . The value of A z will depend 

Z cat Z cat 

on the survey filters, photometric reconstruction technique 
and spectral properites of the survey galaxies. We assume 
a fiducial value of A z = 1 and that the uncertainty on the 
PDF of the population of catastrophic outliers is the same 
as for the slightly scattered galaxies 

cr(P,±J = 6 M (l + z± t ). (23) 



Fig. [3] shows the variation in FoMs with S z and N z , 
the number of tomographic redshift bins, for different com- 
binations of observables, with and without IAs. There is a 
parallel between varying S z and N z , lower redshift error fo- 
cuses the redshift distribution per bin into a tighter shape, 
allowing us to learn from changes with redshift over shorter 
intervals. Somewhat equivalently, increasing the number of 
tomographic bins gives us, in principle, the opportunity to 
detect changes with redshift at higher resolution. 

Increasing 8 Z decreases the constraining power of a sur- 
vey because the loss of redshift precision "smears out" the 
redshift distribution of our sample. In particular galaxies 
leak between tomographic redshift bins meaning we can ex- 
tract less information from tomography. The impact of in- 
creased 5 Z on constraining DE is relatively weak. For ee alone 
the DE FoM drops by only 1/4 from S z = (spectroscopic, 
i.e. perfect redshift precision) to S z = 0.2, already a much 
higher value than the fiducial redshift uncertainties of typ- 
ical Stage-Ill (S z = 0.07) and Stage-IV (S z = 0.05) experi- 
ments. MG FoMs display very similar behaviour. 

The reasons for this behaviour are relatively well under- 
stood. The lensing integral (eqn. [TJ modulates the 3D mat- 
ter power spectrum by the lensing weight functions Wi(x) 
of the redshift bins being correlated. These lensing weight 
functions are broad in redshift and act as kernels smoothing 
the redshift information. This effect limits the usefulness of 
increasingly accurate redshift information and explains why 
the trend with 8 Z is relatively shallow. 

The inclusion of galaxy position information produces 
slightly steeper trends with increased S z because their an- 
gular power spectra integrals are modulated by the redshift 
distribution rii(x) of the bin(s) in question. This is a nar- 
rower function than the lensing weight function, providing 
more scope for improvement in accuracy with better redshift 
information. 

The slightly steeper curves for the MG FoMs may come 
from the larger number of parameters that need to be con- 
strained in this model, which can be done with new types 
of information such as increased redshift resolution. Broadly 
the results for modified gravity are very similar to those for 
dark energy. 

All of the data combinations with varying number 
of tomographic redshift bins, N z , show the characteristic 
plateaux behaviour as N z increases. Without IAs, increas- 
ing the number of tomographic bins becomes much less ef- 
fective after N z ~ 2 — 3. After this number ee is effectively 
flat, while ee + nn and ee + ne + nn flatten out slowly, gaining 
little benefit after N z ~ 10 — 12. This behaviour is expected 
as the lensing kernel and photo-z scatter smear-out redshift 
information and limit the ability of ever-finer tomography 
to improvement our knowledge of redshift behaviour. 

The inclusion of IA terms for the same probes shows 
the expected decrease in constraining power as we have to 
marginalise over IA bias terms. As we have seen previously, 
the MG FoMs take more of a hit due to the IA effects than 
their DE counterparts. 

The decrease in FoM with increasing S z is significantly 
more pronounced once IAs are included, with ee including 
IAs losing 67% of its constraining power for DE over the 
interval probed. Qualitatively we might hope to understand 
this from the fact that the cosmic shear and each IA term, 
II and GI, have a different redshift dependence. Improved 
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Figure 3. Figures showing the variation with photometric redshift error (S z ) [left panels] and number of tomographic bins (N z ) 
[right panels] of a Stage IV survey FoM for dark energy (DE) parameters t»o,tii a [blue lines], and modified gravity (MG) parameters 
Qoi Qo(l+^o)/2 [red lines], as described in 1141 1 and 1151 . relative to a 'baseline' FoM [black cross]. All FoMs contain priors from a Planck- 
like CMB survey. The figures show relative FoMs for an optimistic scenario, in which unccrtainities in the IA model have been excluded, 
[top panels] and a conservative scenario, in which uncertainties in the IA model are marginalized over using a = N z = 5 gridded bias 
model [lower panels]. Three data combinations are considered in each panel, shear-shear correlations 'ee' [dotted lines], 'ee+nn' [dashed] 
and when shear-position cross-correlations 'rie' are included [full lines]. For the galaxy position correlations a = N z = 5 gridded bias 
model is used throughout. 



redshift knowledge allows us to better discriminate between 
the cosmic shear signal and IA contamination. This pro- 
duces stronger constraints on the IA nuisance parameters 
and consequently tighter constraints on DE and MG. 

Similarly, the presence of IAs for increasing N z makes 
the plateau behaviour less pronounced and pushes the level- 
ling off of the FoM back towards N z ~ 6 — 10 for all probes. 
The use of smaller slices in redshift space gains informa- 
tion about the z-distribution of the IA signals and aids our 
constraining power. The effect is curtailed by the same fun- 
damental limit of redshift accuracy in the projected angular 
power spectrum. Again the hit on FoM due to IAs is more 
pronounced for MG than DE but the qualitative behaviour 
of FoM with N, in the presence of IAs is the same for both. 

The effect of changing the fraction of catastrophic out- 
lier redshifts in our survey sample is shown in Fig. 2] As 
we would expect increasing the fraction of completely mis- 
estimated redshifts makes for poorer constraints on both 
MG and DE, for similar reasons to the impact of 5 Z . 

All the lines in Fig.[4]are relatively flat, suggesting that 
the projected nature of the angular power spectra and the 
effect of the lensing weight function continue to limit the im- 
pact of redshift information on overall constraints. No probe 



combination loses more than 30% of its DE or MG FoM as 
f ca t increases from zero to 0.2 which is more than twice the 
typical requirement value of projected stage-IV surveys. 

The general pattern that MG FoMs decrease slightly 
more strongly than their DE equivalents with the addition 
of IAs is repeated here. One interesting feature is that the 
MG constraints for ee without IAs and all the probe com- 
binations with IAs decrease more steeply with increasing 
feat than their DE counterparts. The difference is most pro- 
nounced for ee with IAs where the MG constraints falls by 
twice as much over the f cat range as the DE FoM. 

3.3 Redshift Priors 

As well as describing the change in FoMs with the size of the 
photometric uncertainty, it is important to the quantify the 
effect of our degree of knowledge of that uncertainty on our 
FoMs. In our case we follow the standard procedure of in- 
serting this degree of knowledge as a prior on the parameter 
S z . 

To test this assumption we extend our photometric red- 
shift model- breaking th e redshift range -3 into 30 bins in 
z of width 0.1, following [Ma et al.l (|2006t ). Note that these 
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Figure 4. Figures showing the variation with fraction of catas- 
trophic redshift outliers (f ca t) of a Stage IV survey FoM for dark 
energy (DE) parameters too, fa [blue lines], and modified grav- 
ity (MG) parameters Qo,Qo(l + Ro)/% [red lines], as described 
in J 141 and i!5l , relative to a 'baseline' FoM [black cross]. The 
figures show relative FoMs for an optimistic scenario, in which un- 
certainities in the IA model have been excluded, [left panel] and 
a conservative scenario, in which uncertainties in the IA model 
are marginalized over using a = N z = 5 gridded bias model 
[right panel] . Notation is the same as in Fig. [3] 



bins are used only to parameterise our photo-z model, we 
continue to use the standard 10 bins for shear tomography. 
We allow 8 Z to vary independently in each of these z-bins 
and introduce a new parameter Zbias, the amount by which 
the the photometric distribution of redshifts is offset from 
the true redshift. Zbias is allowed to vary independently in 
each of our 30 z-bins. Each 8 Z has the standard fiducial value 
of 0.05 and each Zbias has fiducial value zero. We have en- 
larged our Fisher Matrix by 60 new free parameters, all of 
which are marginalised over in the results given in Fig. [5] 

We assume that the Gaussian prior assigned to each of 
the 60 free parameters is the same and we allow the prior to 
vary from a width of 0.0001 to 10. Clearly the constraining 
power a given survey achieves decreases as we reduce our 
amount of prior knowledge about the redshift distribution. 

Apart from the expected decrease with wider priors we 
see two distinct regimes- relatively flat FoMs below prior 
0.001 and above 1 with a transition regime in between. It is 
reassuring to note that we cannot achieve arbitrarily large 
FoMs with tighter redshift priors- this is connected to the 
previously noted "smearing" of redshift information in the 
shear angular power spectrum. At large values of the prior 
width we reach a "self-calibration" regime in which the infor- 
mation in t he survey is sufficien t to constrain the unknown 
parameters (|Huterer et al.ll2006h . 

When IAs are ignored, the decrease in constraining 
power is more pronounced. Mostly the lines start higher for 
narrow priors and converge to a similar low FoM for wide 
priors. It is probable that this behaviour is an artefact of 
there being fewer nuisance paramters when IAs are ignored. 
This produces strong constraints in the narrow prior regime, 
while for wide priors redshift uncertainty becomes the dom- 
inant source of error and probe combinations with/without 
IAs are more similar in terms of constraining power. 

What is also clear is that the DE FoM falls off much 
more steeply over the prior range than the MG FoM. This is 



true for all probe combinations, with and without IAs but is 
most spectacular for the lines which include IAs, here there 
is actually a cross-over, with the MG FoM ending higher 
than the DE FoM at prior = 10. For ee, including IAs, we 
observe a a factor of 10 decrease in DE FoM but a factor of 
3 decrease in MG FoM. This trend is repeated when all data 
is used (ee + ne + nn) . We expect the dark energy parameter 
w a to be particularly sensitive to incomplete redshift infor- 
mation because it describes the time evolution of the dark 
energy equation of state. By contrast, the modified grav- 
ity parameters are more robust to poor redshift knowledge 
because they are constant with redshift and the redshift evo- 
lution of the modified gravity functions is fixed by our choice 
of redshift power law s. 

We can relate the required prior on redshift parame- 
ters to the number of spectroscopic quality galaxy redshifts 
needed to calib rate our photometric sample as given by 
iMa et al.l (|2006h : 



28l{l + zf 



(24) 



For ee, ignoring IAs, a prior of 1x10 3 is sufficient to 
recover 90% of the peak FoM for either DE or MG. With a 
photometric scatter of a z = 8 z (l + z) = 0.05(l + z) and a me- 
dian redshift of z m = 0.9/\/2, this translates into 1.34x10* 
spectroscopic galaxies per redshift bin. Including all galaxy 
position information ee + ne + nn and IAs reduces the re- 
quired prior to ~ 1 x 10~ 2 , or of order 100 spectroscopic 
galaxies per bin. This represents a factor of 100 reduction 
in required numbers compared to the standard result. How- 
ever, note that there are other factors which must be taken 
into consideration when determining the number of galaxies 
per redshi ft bin, especially cosmic variance (e.g. see the dis- 
cussion m I Van Waerbeke et al.l 120061; llshak fc Hiratal [iooBI ; 



iBordoloi et alj|20ld : ICunha et a.l.ll201lF 

Fig. [6] shows similar results for the fiducial stage III 
survey. Qualitatively the results are similar to the Euclid 
survey. In this case, with a photometric scatter of a z = 
0.07(1 + 2) and a median redshift of z m = 0.8/\/2, ee without 
IAs requires a prior of ~ 1 x 10 -2 which translates into 
around 100 galaxies per bin. Again ee + ne + nn with IAs 
requires a factor of 100 fewer galaxies, a substantial saving. 
Interestingly, while the MG constraints produce the same 
saving in spectroscopic galaxies required when IAs and all 
probes are included, the absolute values are a factor of 2 
tighter in prior terms (i.e. a factor of 4 more spectroscopic 
galaxies are required). 



4 CONCLUSIONS 



Previous work (|joachimi fc Bridlel2009l ) has shown how con- 
straints from cosmic shear are significantly degraded by the 
inclusion of I A systematics, using a comprehensive grid of 
k- and z-dependent nuisance parameters. The impact of IAs 
can be mitigated by the inclusion of galaxy position corre- 
lations and shear-position cross-correlations using only the 
photometric information gathered as part of the shear sur- 
vey. 

In a related paper, LBKB, we confirmed and extended 
these results showing that using the full ee + ne + nn infor- 
mation with up to 36 nuisance parameters can regain the 
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Figure 5. Figure showing the variation with the priors on the 
bias on median redshift, z^ iaa , of 30 redshift bins spanning the 
survey redshift range due to the effect of photometric redshift un- 
certainties and the photometric redshift statistical dispersion, 8 Z . 
All other survey parameters arc fixed at their fiducial values. Re- 
sults are for a Stage IV survey FoM for dark energy (DE) param- 
eters WQ,w a [blue lines], and modified gravity (MG) parameters 
Qo- Qo(l+-Ro)/2 [red lines], as described in I lf4ll and (15) , relative 
to a 'baseline' FoM [black cross] . As shown in the key in the fig- 
ure, shear-shear correlations, 'ee' , and combined shear and galaxy 
position information including cross-correlations 'ee+ne+nn' are 
shown when uncertainities in the IA model have been included 
and excluded. 



constraining power of a naive ee analysis which ignores IAs 
and that this is true when constraining deviations from GR 
as well as dark energy. In this paper we confirm that these 
results are robust across a wide range of survey parameters, 
encompassing fiducial stage III and stage IV cosmic shear 
surveys. 

Increasing survey area improves both DE and MG 
FoMs, as one expects, whether we ignore or include intrinsic 
alignments. Previous work, which did not include the CMB, 
found a linear variation of FoM with survey area even when 
IAs were included. We have shown, however, that if one 
includes priors from upcoming CMB experiments that this 
dependence is weakened to vary roughly as a power law with 
an exponent of 0.4. The inclusion of IAs slightly weakens the 
relationship between FoM and survey area. 

More realistically, survey time is a constant so increased 
area must be traded off against depth. We find the depen- 
dence on area significantly flatter when IAs are included. 
This is to be expected because IAs dominate over the cosmic 
shear signal at low redshift and we effectively marginalise 
out information from this regime. On increasing the area, 



Figure 6. As in figure[5] but for a Stage III survey specification. 



and thus decreasing the depth to lower redshift, an increas- 
ing proportion of the signal is lost by the marginalisation. 
The effect is more pronounced for a stage-Ill survey than a 
stage- IV survey, which is expected because the depth is very 
shallow indeed if it has to cover half of the sky. 

The quality of photometric redshift information is cru- 
cial to the success of future cosmic shear missions. We show 
that the inclusion of IAs and the full use of galaxy posi- 
tion information requires more accurate redshift information 
than a naive cosmic shear only analysis. GG-only retains 
97% of its constraining power for perfect redshift knowledge 
with a photometric redshift scatter of 5 Z = 0.05, the fidu- 
cial value for a stage-IV survey. By contrast the combined 
ee + ne + nn including IAs requires <5 Z = 0.015 to retain 
the same proportion. Similarly the characteristic plateau be- 
haviour, where WGL constraining power does not improve 
over a tomographic resolution of ~ 4 — 5 bins is pushed 
back on the inclusion of IAs and galaxy position informa- 
tion to ~ 10 — 12 redshift bins. This is understandable as 
the standard cosmic shear redshift dependence, modulated 
by a broad lensing kernel for each redshift bin, has been 
joined by galaxy position correlations which depend on the 
narrower redshift distribution function and the attempt to 
constrain IAs signals, each with their own redshift depen- 
dence. Taken together these place increased importance on 
accurate redshift knowledge. 

We have presented results showing the effect of vary- 
ing our degree of prior knowledge of an extended parame- 
terisation of redshift information, allowing 60 free parame- 
ters to be shared between redshift dispersion and redshift 
bias. Around Azbias = A<5z = 0.1 the probes enter a self- 
calibration regime and the FoMs flatten off. The DE FoMs 
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fall off much more steeply and relatively lower than their MG 
counterparts, highlighting the sensitivity of the DE equation 
of state, particularly w a , to incomplete redshift information. 
By contrast the MG parameters have a fixed redshift depen- 
dence in this work. 

This work assumed a relatively simple model for devi- 
ations from GR. We vary two MG parameters, each with 
the same fixed redshift dependence. In the future we hope 
to extend these results to take into account a MG gravity 
parameterisation which could vary more generally as a func- 
tion of scale and redshift. In contrast our fiducial parameter- 
isation of galaxy bias and I A nuisance parameters is robust, 
containing over 100 nuisance parameters. With a greater un- 
derstanding of the physics of galaxy formation mechanism 
responsible for IAs it may be possible to motivate a more 
targeted, simpler parameterisation. 
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